History effects and pinning regimes in solid vortex matter 
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We propose a phenomenological model that accounts for the history effects observed in ac sus- 
ceptibility measurements in YBa2Cu307 single crystals [Phys. Rev. Lett. 84, 4200 (2000) and 
Phys. Rev. Lett. 86, 504 (2001)]. Central to the model is the assumption that the penetrating 
ac magnetic field modifies the vortex lattice mobility, trapping different robust dynamical states in 
different regions of the sample. We discuss in detail on the response of the superconductor to an ac 
magnetic field when the vortex lattice mobility is not uniform inside the sample. We begin with an 
analytical description for a simple geometry (slab) and then we perform numerical calculations for 
a strip in a transverse magnetic field which include relaxation effects. In calculations, the vortex 
system is assumed to coexist in different pinning regimes. The vortex behavior in the regions where 
the induced current density j has been always below a given threshold (j^ ) is described by an elastic 
Campbell-like regime (or a critical state regime with local high critical current density, j^). When 
the VS is shaken by symmetrical (e.g. sinusoidal) ac fields, the critical current density is modified to 
jf < Tc ati regions where vortices have been forced to oscillate by a current density larger than j% ■ 
Experimentally, an initial state with high critical current density (j^) can be obtained by zero field 
cooling, field cooling (with no applied ac field) or by shaking the vortex lattice with an asymmetrical 
(e.g. sawtooth) field. We compare our calculations with experimental ac susceptibility results in 
YBa2Cua07 single crystals. 

PACS numbers: 74.60. Ge, 74.60. Jg 



I. INTRODUCTION 

The competition between vortex-vortex interactions 
and randomness determine the electromagnetic_response 
of type-II superconductors in the mixed state.EJ The in- 
teraction between vortices favors an ordered vortex lat- 
tice (VL) that opposes the disorder arising from random 
pinning centers and temperature. The combined effect of 
these interactions leads to a great variety of collective be- 
haviors, making the VL a paradigmatic system to study 
an elastic object in a disordered environment .el 

The rich static phases of vortex matter in high T c su- 
percondurjtors have drawn a great deal of attention in 
the pastjjtrcl but recently, much effort has been directed 
towards understanding vortex lattice dynamics both in 
low and high T c materials. Though details are not com- 
pletely understood, consensus has been reached that the 
nucleation and annihilation of defects (e.g., dislocations) 
in the VL might 1 Blay 1 a ma jQE-xoleJp its response to an 
applied forcelSQffiBBBiyye'N Reproducible states 
with distinct mobilities have been observed, and it has 
been argued that this reflects distinct degrees of topolog- 
ical order of the VL: a defective VL would interact more 
efficiently with pinning centers than an ordered one be- 
cause of a reduction in its effective shear modulus. In 
particular, several experiments have revealed that, in cer- 
tain conditions, an ac maso^ju^fifiH^paji-ac current can 
assist the VL in orderinglBBBMBO 

At first, two basic mechanisms have been invoked to 
account for the dynamic ordering of the VL in low T c su- 
perconductors: a dynamic transitio»,pccy|rring at current 
densities above a certain thresholdalla'B and an equili- 



bration process where an ac magnetic field (or current) 
perturbation assists the VL in going from a supercooled 
disordered state into the more ordered equilibrium state. 
However, data have been reported that can not be ac- 
counted for within either of these two frameworks. We 
refer to the recent findings in the vortex lattice in twinned 
YBa 2 Cu 3 7 (YBCO) single crystals.BB It was found 
that when vortices are shaken by a temporarily symmet- 
ric ac magnetic field (e.g., sinusoidal), they are driven 
into a mobile state, but if the shaking field is temporarily 
asymmetric (e.g., sawtooth), the lattice is trapped into 
a more pinned structure.B As suggested in Ref. [D| the 
possibility of switching to a high mobility VL state by 
applying a symmetrical (sinusoidal, square, triangular, 
etc.) ac magnetic field or to a low mobility VL state by 
applying a temporarily asymmetric magnetic field wave- 
form (sawtooth, etc.) is inconsistent with a naive equili- 
bration process. In addition, the induced current density 
that participates in ordering or disordering the VL can 
be tuned to be comparable by adjusting the magnetic 
field amplitude and frequency. Disordering observed for 
higher field sweep rates than those that are used in sym- 
metric oscillating fields that order the vortex lattice is 
a clear indication that the higher mobility of an ordered 
VL is not attained-through a dynamical transition at high 
current densities .B If this were the case, the VL should 
also order with the asymmetric ao-field that is expected 
to induce a higher current density.E3 This demonstrates, 
at least for YBCO, that the high or low VL mobility 
results from a collective dynamical process of a different 
kind that we believe is intrinsic to the oscillatory motion. 

The attained degree of mobility is found to persist af- 
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ter interrupting the ac excitation for long periods of time 
(exceeding 1 h). As the ac field penetration is depen- 
dent on its amplitude, Hq, and temperature, T, different 
regions of the sample may remain in different pinning 
regimes depending an the detailed thermomagnetic his- 
tory of the samplco The dynamical rearrangement of 
the VL occurs in the outer zone of the sample where the 
penetration of the ac magnetic field induces ac currents 
that force vortices to move. An increment of Ho or T will 
result in an increase in the penetration depth of the ac 
field. As the ac field further penetrates into the sample, 
the VL mobility changes. In Ref. [lij we argued that, as 
long as the current density does not exceed the critical 
one (so that a small fraction of the vortices moves out of 
their pinning sites) the VL will remain in the low mobil- 
ity state. In a bulk-pinning-dominated regime, this will 
happen in the inner nonpenetrated portion of the sample. 
As T or Hq increases this region will shrink and, eventu- 
ally, disappear. The robustness of the states involved and 
their spatial distribution in the sample are the main in- 
gredients leading to the observed memory effects in Ref. 

E 

In this paper, we present analytical and numerical cal- 
culations for the field and current distributions, as well 
as for the ac susceptibility Xac in samples with inhomo- 
geneous pinning. If history effects and flux creep can be 
neglected, the current and flux densities inside the super- 
conductor a#« often well described by the Bean critical 
state modelB (CS) which assumes a critical current den- 
sity that is magnetic field independent and a null lower 
critical field, H c \ (H c \ = 0, B = hqH). Exact analyt- 
ical solutions, [have b een obtained for the perpendicular 
geometry 2 2 H 2 'rr 4 lr 5 rr 6 l which is more appropriate for our 
experimental arrangement. For more general situations, 
as in our particular case, a time-dependent theory is re- 
quired. If it is assumed that the critical current density 
does not vary over the sample thickness, this dynamic 
theory can be formulated in terms of a one-di men sional 
integral equation as discussed by E. H. Brandt El cB With 
the use of an appropriate constitutive relation E(J) it is 
possible to compute the flux motion (and the magneti- 
zation) in strips, disks, etc. with homogeneous or in- 
homogeneous j c when the applied field is cycled. From 
the magnetization curves, we calculate the susceptibility. 
These results are compared with measured ac suscepti- 
bility data obtained in a twinned YBCO crystal following 
a specific protocol to control the spatial pattern of the 
dynamical behavior of the VL in the sample. 

As an example, the VL can be prepared by zero ac 
field cooling (ZF ac C), cooling the sample below the crit- 
ical temperature before applying the masaetic ac field. 
This traps the VL in a low mobility statetil with a high 
critical current density, ~ j> . If a sinusoidal ac field 
is then applied, the VL will become more mobile in the 
penetrated outer zone of the sample and it will be charac- 
terized by a lower critical current density, ~ . Vortices 
in the inner region of the sample, where the induced cur- 
rent density j has been always below j> , will remain in a 



more pinned structure in which vortices oscillate confined 
in their pinning potential wells. 

The critical current patterning so defined may be ob- 
served, for example, by removing the ac field and apply- 
ing an increasing ac probe to measure the ac susceptibil- 
ity Xac- This is the basic protocol that we have developed 
in calculations and in our ac susceptibility measurements 
to show evidence of the different local degrees of VL mo- 
bility. Analytical calculations that qualitatively intro- 
duce our model are presented in Section II A. In Section 
|IIB| we present our numerical calculations of \ac that 
take into account a more adequate transverse geometry 
and also include the effect of flux creep. In calculations, 
we consider two possible values for the critical current 
density though it should be kept in mind that it is ob- 
served experimentally that j c depends on the number of 
cycles of the "shaking" ac field (see Ref. |l5|) . Experimen- 
tal results arc described in Section III and conclusions are 



drawn in Section [V 



II. THEORY 

For symmetrical shaking fields, we can argue that the 
VL in the inner portion of the sample will remain in a 
low mobility state as long as the inner current density 
does not exceed the critical current density. In this case, 
a Campbell-type regime may describe the elastic pinning 
of vq rt i ccys tha t-arc confined to oscillate inside a potential 
welLBBElM 

The main criterion to discriminate between a Camp- 
bell and a critical state regime arises from comparing the 
scale of the vortex displacement, Ax, produced by the in- 
duced currents, and the size of the pinning potential well, 
r p . The application of an ac magnetic field results in a 
periodic compression of the VL, with a maximum effec- 
tive vortex displacement at the sample's boundary. This 
displacement can be easily estimated from the variation 

of the lattice parameter a = J -j^bJ~ t^ a t occurs after 
applying a perturbation in the dc field, Bd c - For a sam- 
ple with no pinning and with [i§H = B (i.e., neglecting 
H c i), the vortex displacement at the sample perimeter 
arising from an ac field of amplitude Bq (Bq << Bd c ) 
is Ax « i-J^-i? where R is the sample radius. If Ax 
is smaller than r p the elastic pinning regime usually de- 
scribes satisfactorily the penetration of magnetic flux. As 
the ac field amplitude increases, the region in the sample 
where the elastic regime applies naturally shrinks. We 
shall come back to this criterion further on, and momen- 
tarily we shall address the problem of a superconductor 
with nonuniform critical current density. 



A. Statics. Analytical calculations 

We start by calculating the ac susceptibility of an 
infinite slab (|x| < A) in a longitudinal magnetic field 
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FIG. 1: (a) Infinite slab for z and y used in the estimation 
of the ac magnetic susceptibility. The field is applied in the z 
direction. The critical current density is inhomogeneous, with 



for \x\ < Xd (dark gray) and j c 



for 



\x\ > x d 



(light gray), (b) Field profiles inside the sample for — > oo. 



for which the inner critical current density, j % c nt = j> 
(\x\ < x c < A) is larger than the outer one, j% xt = j< 
[x c < \x\ < A) [see Fig. 0(a)]. This first calculation will 
provide a qualitative description of the observed ac field 
amplitude dependence of the ac susceptibility. We shall 
not consider at present that the ac field modifies the 
local critical current density, and we will suppose a given 
fixed spatial distribution of j c . For illustrative purposes 
we also take j> — > oo so that the outer region of the 
sample Xd < \x\ < A is in the CS whereas there is perfect 
shielding (Campbell penetration length Ac -C X4) for 
\x\ < Xd- The calculation of \ ac in this simplified 
situation is straightforward and fallows from known 
results for the critical state in a slato if one notes that: 
i) for ac field amplitudes below Hd 



j<{A - X d ), Xac 

is the well known CS susceptibility for an infinite slab 
with j c = jf; ii) for fields higher than Hd, the sample 
magnetization is given by the sum of the magnetization 
due to the outer region of total width 2(A — Xd) in 
the CS and that due to the inner region of width 2xd 
with perfectly diamagnetic susceptibility x' — *x" = — 1 
[see Fig. [ljb)]. For an applied ac magnetic field, 
H ac = Hq cos(wi), the susceptibility is 
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As an example, we plot in Fig. ||(a) the calculated 
imaginary component x" as a function of the normal- 
ized ac magnetic field Hq/ H p for different values of Xd/A. 
These results describe the case of a fixed boundary sep- 
arating high and low mobility regions of the VL (in each 
curve the boundary is at a different position). We will 
address now the more realistic situation where the ac field 
modifies the position of the boundary. 

As was discussed above, the elastic pinning regime is 
valid for ac fields of low amplitude, where vortex displace- 
ments at the sample boundary do not exceed the size of 
the pinning potential well. As the ac field amplitude is 
increased, the region in the sample where this condition 
is valid will be gradually reduced. As a consequence, the 
curves plotted in Fig. @(a) for fixed Xd can not be valid 
for an arbitrary high ac field amplitude. A more realistic 
approach has to consider that Xd will gradually reduce 
for sufficiently strong fields. This reduction should be 
monotonic in the ac field amplitude. 

It is useful to define a "minimal penetration suscepti- 
bility" that represents the ac susceptibility of a sample 
that initially contains a homogeneous low mobility VL 
and is subjected to an increasing ac magnetic field start- 
ing from zero. These curves are represented in Fig. |^(a) 
by dotted lines and will be discussed in the following 
paragraphs. 









0.2 










b) ^^^^ 


0.1 






0.0 







A-x d - 



FIG. 2: Dissipative component of the ac susceptibility, x" 1 
following Eq. | and [| (a) X " vs. for x d /A = 0, 0.2, 0.4, 
0.6, 0.8 (solid thick- lines) . The arrow indicates increasing 
Xd/A. Xd = corresponds to the usual critical state. Dotted 
lines represent the curves of minimal penetration defined in 
the text for hr = 0.2, 0.4 and 0.6H P . (b) Susceptibility for 
Ht = OAHp and initial Xd = x d /A — 0.4 (solid thick-line). 
The arrows indicate the direction of the magnetic field varia- 
tion. 

The application (in the elastic regime) of a small per- 
turbation field b s (ac or dc) induces currents that exert 
a pressure on the vortices at the sample boundary that 
propagates inside the sample due to the VL elasticity. 
In a wide slab, the perturbation decays exponentially, 
with a chap^ftfijEifttic length Ac, the Campbell penetra- 
tion depth.EE 
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We suppose first an increasing ac field amplitude start- 
ing from H = 0. As H is increased above a given 
threshold, Ht, vortices in the outer part of the sample 
are driven out of the Campbell state as they are forced 
to move out of their pinning sites. The subsequent os- 
cillatory motion of the vortices produces a reordering of 
the VL that results in a local reduction of j c to jf . Here, 
Kt is related to the depth of the pinning potential well 
in which vortices oscillate in the Campbell regime. 

We can estimate Xd as a function of H for a semi- 
infinite sample (x > 0). A perturbation b s penetrates 
the sample as AB(x) = b s exp(— x/X ac ), where the ac 
penetration depth is given iii^X 2 ^ = A 2 + A c with A the 
London penetration depth£3'E3 In the case we are inter- 
ested in Ac = (cn/aL) 1 / 2 , where cu is the compress- 
ibility of the VL and is the Labusch parameter (in 
general for weak pinning Ac ^ A and X ac ~ Ac). 

To determine whether the elastic limit is applicable or 
not, we estimate the vortex displacement at the sample 
surface as 



Ax 



AB(x) 
B dc 



dx 



(5) 



As discussed above, for displacements Ax ~ r p , or 
larger, the elastic regime is not expected to be applica- 
ble. We find from Eq. [5] that this occurs for an applied 
perturbation b s = br ~ fp\/ a L^o, where we have used 
A ac ~ Ac and cu ~ Bj c /fi (kr = Hoh T )- 4f the VL is 



not strongly perturbed then 



'Bdc £3 leading to 
br oc yfj>. 

In the semi-infinite sample, it is evident that this rela- 
tion establishes the applicability of the elastic Campbell 
regime at any point inside the sample. Then, we can as- 
sume that, if the VL is in a low mobility state (J c = j^) in 
the whole sample, and the ac field amplitude is increased 
from zero, Xd will reduce monotonically from Xd = A to a 
value given by AB(xd) = br oc \[J> ' . In this picture, the 
ac susceptibility of a slab is given by Eqs. (l)-(4) with 
Xd not a constant anymore but depending on the applied 
ac field amplitude, 



H <h T - 
H > h T : x d 



Xd 

= A 



A 

Ho 



Fig. ||(a) shows in dotted lines the "minimal penetra- 
tion curve" as defined above for the imaginary compo- 
nent x" as a function of the normalized ac magnetic field 
H / H p for hr = 0.2, 0.4 and 0.6 H p . In Fig. |(b), we rep- 
resent the ac susceptibility of a slab predicted by this sim- 
plified picture where the initial state corresponds to a sit- 
uation in which the vortices in the outer region, |x| > Xd, 
have a higher mobility (j c = j',5). In this example, we 
choose hr — QAH p and we show x" f° r different initial 
values of Xd (x d )- The curve represented with a thick 
line corresponds to x\ = 0.4A For H < Hd the internal 
region of the sample provides no contribution to x" 1 and 



the susceptibility corresponds to a slab in critical state. 
For H d < H < H d + h T (in this example H d + h T = H p ) 
Xac moves away from the CS susceptibility as the internal 
region contributes with a higher screening capability and 
the dissipation decreases. However, as AB(xd) < br, Xd 
remains unchanged and equal to its initial value, x l d . For 
Ho > Hd + Ht{= H p ) the boundary between low and 
high VL mobility progressively penetrates into the slab 
and Xac follows the minimal penetration curve. Finally, 
for Ho = H p + Ht{= lAHp) the sample is fully pene- 
trated and the CS regime is recovered. The convenience 
of this calculation will become clear when we describe in 
Section IV the protocols followed in experiments. 



B. Numerical calculations 

1. Computational method 

The original Bean critical state model has been— ep, 
tended to different-geometries, including long stripscSE^I 
and circular diskst3'Ot3 in a perpendicular magnetic 
field. These geometries are more appropriate for our ex- 
perimental situation. The results obtained analytically 
can be reproduced by an algorithm developed by E. H. 
Brandt EIlEj The great advantage of the numerical ap- 
proach is that it allows to analyze static and dynamic 
problems in samples where the critical current density 
has a non-trivial dependence on position. The supercon- 
ductor is modeled by a constitutive law E(j) that has 
the form E{j) = Eo(j-) p corresponding to a logarithmic 
activation energy U(j) ~ hi(j c /j) or to U(J) - {j c /j) a 
with a <C f for a superconductor with critical current 
density j c . 

Following Brandt, we have considered a strip localized 
in \x\ < d/2 < A with \y\ < A and \z\ > A in an 
applied field H a x and we define the sheet current density 
3(y, z) — J_d 2 /2^( x ' V' z )dx, which in this case is J(y, z) — 
—J(y)z (the extension of these calculations for disks is 
straightforward and is described in Ref. |28| ). 

The magnetic field generated by these currents is given 
by Ampere's law. For the strip, the total field perpen- 
dicular to its surface evaluated at x — has the form 



H x (0, y) = H(y) = ±- [ ^du + H a 



(0) 



As the current flowing in the strip is induced by the 
external magnetic field, the symmetry J(y) — —J(—y) 
can be used in Eq. ^]. Considering the electric field 
E(y, t) — <f>(y, t)/L, induced by the time dependent mag- 
netic flux 4>{y) = Ho Jo H(u)du, one gets an integro- 
differential equation, for the sheet current. This equation 
has been discretizedEZrEj 
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N 



■/,•(*) 



2iru 3 H a {t) 



(7) 



where r = p^adj (2tt p) , p = E(j)/j and iV is the number 
of points (0 < iij < 1) of the spatial grid, = u(^) 
with ii = (i - \)/N. In Eq. the strip half width, A, is 
the unit of length and is given by (see Ref. gTJ) 



A"; 



Sin 

iV 



Ui+Uj 



i ^ 



(8) 



~i In w j 
iV 111 4,-irNui 



where Wi = u' i is the weight function. 

Having obtained the sheet current as a function of time 
by solving Eqs. @-(@), we evaluate the magnetic mo- 
ment m of the sample that points in the x direction as 



§ / j x r d 3 x 



which for the strip results in 



m = -xML, M = f_ A yJ(y)dy 



(9) 



(10) 



If the applied field is sinusoidal, H a (t) — H ac (t) = 
Hocos(uit) where u> — 2ir/T, it is possible to define the 
ac susceptibility coefficients as 



creep and p > 50 would require long simulation times to 
detect relaxation effects. Note that, in general, only the 
combination E c / J? is relevant (Eq. |l3|), and the voltage 
criterion, E c , and J c are not independent. While in the 
critical state limit E c is irrelevant, in the ohmic regime 
E c /J c = P is the resistivity. Within the dynamical ap- 
proach that we describe, J c has a different meaning. As 
the current is time dependent (it may depend on the rate 
of variation of the magnetic field) , J c is related to the de- 
gree of vortex mobility, where low values of J c represent 
highly mobile vortex structures. 

The space variable was discretized with N — 50 
nonequidistant points (results with N = 30 showed no 
significant difference). The tabulation was performed so 
that the weight function w(x) = u'(x) vanished at the 
integration boundary to remove the infinity in the inte- 
grand occurring at u — > 1 at t = (see Ref. |27| ) 



u(x) 



w ( x ) = -ji)--* 1 ) 

For the ac field H ac (t) = HQCOs(u>t) measured in units 
of j c d, lo = 2E c /(po j c dA) which, in units of A = 2E C = 
j c d = 1 is u> = 1/jUQ. The ac susceptibility was calculated 
after two complete ac cycles that are necessary to reach a 
stationary state in the calculation of the hysteresis loop 
M[H(t)]. 



2. Current and field profiles, ac susceptibility 



w r 

Xn = ^T M(t) cosintut) dt (11) 
^n = -FT M(t)sm{ncjt)dt (12) 

To obtain the susceptibility we solved Eqs. (@)-(||) by 
a fourth order Runge^Kutta method with automatic con- 
trol of the step size.E3 We chose the constitutive relation 
suggested above, 

E(J,= Et (-Q P i ,13) 

with J = | J|. 

In the Bean critical state model E(J) is strongly non- 
linear and has an abrupt increment at J = J c . For the 
constitutive relation of Eq. |l3| this occurs for p — > oo. 
For 1 <C p < oo, it leads to a nonlinear relaxation (ap- 
proximately logarithmic) of M(t) caused by the nonlinear 
magnetic flux diffusion. As p decreases the relaxation 
becomes stronger and in the limit p = 1 the behavior 
is ohmic, with linear magnetic flux diffusion and expo- 
nential relaxation of M(t). In our simulations we chose 
p = 20, as p < 10 would blur the flux front due to flux 



In our calculations, the critical sheet current J c at each 
point j/i inside the strip depends on the magnetic history 
of the sample. As was discussed above, we assume that 
J c can take two different values, J> or Jf (now both fi- 
nite) depending on the sheet current J(yt) that flowed in 
t/j in the previous iteration. For a given Hq, if J(t/j) > 
then J c {y{) is replaced by Jf- (which has been arbitrarily 
chosen as unity) and the calculation is restarted. Clearly, 
the magnetic history of the sample fixes the initial dis- 
tribution of J c (y)- In Figs. Hand[§ we show the result of 
our calculations starting with zero ac field and different 
initial conditions. In Fig. [| we show the field (upper 
panel) and the current (lower panel) distributions for a 
strip with the following initial critical sheet current dis- 
tributions: J c —Jf (left), J c —J^ (center), and J c =J< 
for \y/A\ > 0.75 and J C =J> for \y/A\ < 0.75 (right pan- 
els). 

In the protocol for the calculation of ac susceptibility 
(Fig. |^) we sweep the ac field from zero. The minimal 
penetration curve is calculated as follows. We begin with 
J c = J> for the whole strip (low mobility initial state, 
center panels in Fig. ^). The field is increased from zero 
and we compute the field and current profiles. If the 
sheet current exceeds Jf at any point of the strip, then 
J c is changed to J^ at that point and the calculations 
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FIG. 3: Field profiles (top) and current profiles (bottom) obtained numerically for a field H cos cut increasing from zero and 
evaluated at Ho, for different Ho < 1. Left: homogeneous sample, = Jf = J c . Center: penetration for the minimal 
penetration curve = 2J,5. Right: profiles obtained increasing Ho with initial x r d = 0.75A Calculations are for a thin 
superconducting strip. See the text. 



are restarted. This generates by construction the small- 
est outer region with J c = J< for each given ac amplitude 
Hq (so that J < J> in the whole sample). The field am- 
plitude Ho is increased again and the procedure repeated. 
Note that the region with J c — J> shrinks towards the 
center of the strip as Hq is increased. In the calculation 
of Xac the field has to be varied in a complete cycle, but 
it is clear that the critical current density distribution is 
fixed after the first quarter of cycle. 

Now we calculate the rest of the curves that result for 
a strip with an initial distribution of J c as in the right 
panels of Fig. ^. To obtain the initial state, we apply a 
magnetic field of a given amplitude and we define the ini- 
tial distribution of J c following the procedure described 
above. We then calculate the ac susceptibility increas- 
ing Hq from zero and check for the condition to replace 
J> by Jf. This occurs for field amplitudes larger than 
the one that defined the initial state. Note that from this 
field amplitude on, the calculated susceptibility coincides 
with the minimal penetration curve. 

It is worth noting that in a strictly 2D model, it is 
not possible to have an abrupt change in J c , because this 
causes a local divergence in J and vortices would always 
move. To solve this inconvenience, in our calculations, 
the critical current density changes smoothly from J> 
to J< over distances Q-f-ithe order of d (where the 2D 
treatment breaks downed). 



In Fig. | (left panels) we took J> = 2J<. In Fig. 
^ (right panels) we compare the minimal penetration 
curves for different values of J> (1, 1.5 and 2). As we 
shall see in the next section, the protocol followed in the 
numerical calculations is particularly convenient to di- 
rectly compare the numerical results with the experimen- 
tal ones. In the experiments that are described below, we 
are able to prepare, in a highly reproducible way, an ini- 
tial state with low mobility in the whole sample (that 
is, described in calculations with J>) or an initial state 
with high mobility at the outer region of the sample (a 
situation that corresponds to J< in the outer region and 
J> in the inner region). 



III. EXPERIMENT 

We measured the ac susceptibility of a twinned YBCO 
single crystaEa using a standard mutual inductance tech- 
nique. The characteristics of the crystal are: dimensions 
0.56 x 0.6 x 0.02 mm 3 , T c = 92 K and ST C = 0.3 K de- 
termined by ac susceptibility (h ac =l Oe) at zero dc field. 
In our measurements the ac field is applied parallel to 
the c axis of the sample. The dc magnetic field H^c = 3 
kOe was applied at 20° out of the twin boundaries to 
avoid the Bose-glass phase. Further details can be found 
in Ref. 0. 
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H /J *d H n /J *d 

c c 

FIG. 4: Numerical calculation of the ac susceptibility of a strip, \" (top) and \' (bottom). On the left, we show results 
increasing the field for different initial values of Xd, x d /A = 0, 0.05, 0.25, 0.5 and 1. = 2,7,5. On the right, we show the 
minimal penetration curves for J> = 1, 1.5 and 2 Jf. They are obtained by taking x d equal to zero. See the text. 



15). This field 



The generally unknown dependence of the critical cur- 
rent on temperature makes \ac measurements vs. ac 
field amplitude at constant temperature the best choice 
for comparison with numerical results. The experimental 
procedure, based on our previous results, is as follows: 

• We set the VL in a low mobilty state by applying 
10 4 cycles of an ac magnetic field with a sawtooth 
waveform (10 kHz, 8 Oe, see Ref. 
completely penetrates the sample. 

• We define the initial state by applying 10 5 cycles 
of a sinusoidal ac field of a given amplitude at 10 
kHz. 

• This ac field is turned off and we measure ac sus- 
ceptibility increasing the amplitude of the ac field 
probe from zero, with a fixed measuring frequency 
of 10 kHz. 

Results are shown in Fig. [s] for a fixed temperature 
and Hdc = 3 kOe. Different initial states were prepared 



applying 10 5 cycles of a sinusoidal ac field of amplitudes 
0, 1.6, 2.4, 3.2, 4, 4.8 and 8 Oe. 

The similarities with numerical simulations for J> = 
2J< shown in Fig. |] are evident. The marked change 
in the response at the ac field amplitude that defines 
the initial state supports the scenario of vortices being 
shaken or not out of the pinning centers on either side 
of a well-defined boundary. The effective critical current 
density of the more mobile state can be estimated from 
the position of the dissipation peak. Comparing Figs. |] 
and | we have j< ~ ^ eafe /(0.75d). With d = 30 Aim and 
Hl eak = 3.5 Oe, j< ~ 1250 A/cm 2 , and thus j> ~ 2500 
A/cm 2 (at 85.2 K). 

We also observe that the inductive component does not 
approach -1 as Hq is reduced to zero. This is probably a 
consequence of the finite penetration depth in the Camp- 
bell regime £3c3 We can determine some of the relevant 
parameters in this regime, estimating from the satura- 
tion value (x'(0) ~ -0-9), the value for Ac ~ A ac ~ 10/zm 
(^> A ~ 0.6/zm, Ref. ^lj) and the size of the potential well 
^ach-T I Hdc with hr ~ 1-5 Oe and 



50 A (Eq. | ■ 
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FIG. 5: Experimental results, x' an d x" f° r a YBCO sin- 
gle crystal as a function of the ac field amplitude, Ho- The 
measurements were performed increasing Ho from zero. The 
different curves correspond to different initial conditions. Be- 
fore the measurement starts, a low mobility VL is prepared 
in the whole sample and an ac field of amplitude Hq and fre- 
quency / = 10 kHz is applied to define the initial state. Hq = 
0, 1.6, 2.4, 3.2, 4, 4.8, and 8 Oe. 

Hdc = 3 kOe). The value for hx was estimated as the 
minimum ac amplitude capable of modifying the mea- 
sured susceptibility (starting from a low mobility VL). 

We note that from these estimates r p turns out to be of 
the order of the coherence length £ ~ 60 A at our working 
temperature (£ - f (l - T/Tc)- 1/2 with £ - 20 A, Ref. 
fll|). The Labusch parameter is a L ~ ~ ~ 

10 N/m 2 ^ (100 dyn/cm 2 ^). 

3. Frequency dependence of the penetration depth 

Within the dynamical model described above, an in- 
crease in the frequency of the ac field induces a higher 
ac current (on the average) with a smaller penetration 
depth. One can gain some insight by considering a sta- 
tionary state where a constant field ramp H a is applied. 
In the stationary state the current induced by the field 
ramp does not change with time and the electric field 
is obtained as V x E = — B a , where B a = /ioH a , im- 
plying E ~ B a . Inverting Eq. || we get J ~ (B a ) 1/p , 
which shows how J grows as B a is increased. Taking 
J=constant in, for example, Eq. we arrive at the same 
result. 

The reduction in penetration due to an increment in 
the frequency can be deduced using an interesting scaling 
property of Eq. [7] (Ref. . Explicitly, if the unit of time 
is changed by a factor k and time is expressed as t = t/k, 
then the new functions J(t) = J(t) k 1 ^ -1 and B a (t) = 



B a (t) k 1 ^- 1 satisfy the same equations. The resulting 
magnetic field scales with the same factor. If we consider 
a periodic magnetic field, H ac (t) — Hosinut, this scaling 
property implies that if the frequency is increased by a 
factor 10 and the ac field amplitude is also increased but 
by a factor 10 < - 1 / p ~ 1 ' ) , the field penetration depth inside 
the sample will remain unchanged.^ 2 ] This leads to the 
immediate conclusion that if the frequency is increased 
but the amplitude of the field is maintained constant, 
then the penetration depth will be smaller. 

This behavior is observed experimentally. In Fig. ^ we 
compare the imaginary part of the susceptibility mea- 
sured at 10 kHz after applying 10 5 cycles of a sinusoidal 
ac field of a given amplitude at 100 kHz [Fig. ^(a)] and 
at 10 kHz [Fig. ^(b)]. As the amplitude of the mea- 
suring field increases, the susceptibility senses the VL 
mobility in the regions closer to the center of the sample. 
As the change in the response towards the low dissipa- 
tion branch is first observed in Fig. ^(a), it is clear that 
the penetration at 100 kHz is smaller than at 10 kHz. 
The reason for this is that, in the former case, the high- 
mobility outer-region is necessarily smaller as implied by 
the susceptibility measurements (note that the suscepti- 
bility is measured at 10 kHz). From this difference, and 
applying the scaling law discussed above, we can esti- 
mate the exponent p in Eq. [ll^. As an example, the 
penetration depth of 3.2 Oe at 100 kHz is approximately 
equal to the penetration of 2.7 Oe at 10 kHz. Assuming 
a linear dependence of the penetration depth with the ac 
field amplitude, we obtain p ~ 15, which is close to the 
value we used in our simulations. 



4- Temporarily symmetric and asymmetric ac magnetic 
fields 

We have suggested in Ref. [l5] that a change in the mo- 
bility of the VL may be related to a change in its degree 
of order, and that these changes are intrinsic to thp-pscil- 
latory dynamics. Molecular dynamics calculations^ sup- 
port the scenario where defects are healed by a temporar- 
ily symmetric oscillation and created by a temporarily 
asymmetric oscillation. These simulations also show that 
the mean velocity of the VL increases as the number of 
defects decrease. 

In the previous section we showed that it is possible to 
increase vortex mobility (or to order the VL) by apply- 
ing a temporarily symmetric ac field, and we proved that 
this reordering occurs in the outer penetrated region of 
the sample, i.e. in the regions where induced currents 
have densities above the critical current density. We also 
know that a temporarily asymmetric field produces the 
opposite effectj-jt leads the VL into a low mobility (dis- 
ordered) stateJlj Then, the following questions arise: i) 
Does the penetration depth of the disordered region have 
a similar dependence on ac amplitude as the penetration 
depth of the ordered region for symmetrical ac magnetic 
fields? ii) Does the amplitude of the asymmetric field 
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FIG. 6: Experimental results analogous to those shown in 
Fig. ^. To define the initial state, we applied an ac field 
with the amplitudes used in Fig. ^ but with frequencies 10 
times higher (/ = 100 kHz) [(a)]. As a comparison we have 
included results for x" from Fig. |^ in panel (b). It is clear 
that in the latter case the field that defines the initial state 
(10 kHz) penetrates deeper into the sample. This is because 
the boundary between low and high critical current densities 
is found at a higher Hq when the initial state is defined with 
10 kHz as noted by the vertical lines. In other words, the 
vertical lines compare the corresponding ac field amplitudes 
for which the measuring field has penetrated inside the sample 
a length similar to that penetrated by the field that defines 
the initial state. 



have to be above a certain threshold to produce any ef- 
fect? 

Related to this last question we found that, for the ap- 
plication of a symmetrical field, the minimum amplitude 
that produces some reordering is ~ 1.5 Oe, and this value 
may be related to the depinning of vortices out of their 
pinning sites. As we shall see, however, the threshold 
amplitude for the asymmetric field to produce some dis- 
ordering is twice as high, thus the process that prevents 
the VL from switching to a more pinned state should 
have a different origin. 

The experimental protocol for the asymmetric ac fields 
is very similar to the one described above for symmetric 
fields, but was modified so that the mobility of the initial 
state is low. The following procedure is used: 

• We generate a high mobility VL in the whole sam- 
ple by exposing it to 10 4 cycles of an oscillating 
magnetic field with a sawtooth waveform (8 Oe 
amplitude and 10 kHz) followed by 10 5 cycles of 
a sinusoidal ac field of the same amplitude and fre- 
quency. This leads to a high mobility VL in the 
whole sample. 

• We define the initial state by applying 10 5 asym- 
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FIG. 7: x' an d x" m increasing field for different initial con- 
ditions. The mobility of the VL is reduced when asymmetric 
ac fields are applied. Before the measurement starts, a high 
mobility VL is prepared in the whole sample and a sawtooth 
field of amplitude Hq and frequency / = 10 kHz is applied to 
define the initial state. Hq — 0, 1.6, 4, and 8 Oe. 



metric ac magnetic field cycles (i.e., sawtooth) of 
a given amplitude, Hq (10 kHz). 

• We measure ac susceptibility with a sinusoidal ac 
probe as a function of amplitude, starting from am- 
plitude zero (measuring frequency: 10 kHz). 

Results are plotted in Fig. fflfor the same temperature 
and field conditions of Figs. ^] and [| We investigated 
different initial states prepared by the application of a 
sawtooth waveform ac field with amplitude, Hq, between 
and 8 Oe (we show the results for Hq = 0, 3.2, 4, and 
8 Oe). The first observations are that there seems not 
to be a clear boundary separating two regions with dif- 
ferent vortex mobilities as in Fig. g, and that changes 
in susceptibility are important only for Hq ~ 3 Oe or 
higher. With asymmetric ac fields, the curves for the 
different Hq resemble the minimal penetration curves of 
Fig. | for different J l c nt s. The fact that there is no clear 
boundary between two regions with associated different 
mobilities of the VL might be explained by the differ- 
ence in the sweep rate of the sawtooth waveform that 
implies a different penetration depth in the ramp up and 
down of the field. The measuring ac field can also con- 
tribute to blur the boundary as it tends to counteract 
the effect of the sawtooth ac field. On the other hand, 
we cannot obviate the fact that the shape of the suscep- 
tibility does not favor its observation. Finally, another 
possible explanation for this behavior is related to the 
mechanism that leads to the decrease of the VL mobility. 
In Ref. [n| we proposed that this decrease is related to 
an increase in the density of topological defects due to 
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plastic distortions in the VL that are produced by the 
temporarily asymmetric induced currents (cf. Fig. [| it 
is implicitly shown that higher field rates induce higher 
currents). The observation of a minimum amplitude ~ 3 
Oe required to produce a strong effect in the mobility 
of a moving lattice supports this assumption. For plastic 
distortions to occur, we expect that vortices should move 
around one lattice parameter. For a sample in the critical 
state, the displacement at the sample boundary caused 
by a field that penetrates a distance d is Ax w jffi-d 
(for d smaller than the sample radius). For a sample of 
radius 250 /jm in a 3 kOe dc field and 3 Oe ac field (6 
Oe peak to peak, and 5 Oe of full penetration) the total 
vortex displacement at the sample periphery is ~ 0.1 /im, 
which is of the order of the lattice parameter a = 0.09 
/im. 

We note that although an important portion of the VL 
is moving by the application of 3 Oe ac field, only the 
portion near the sample boundary is initially distorted. 
This is not what occurs during VL ordering, where all the 
vortices that move by the action of a symmetrical ac field 
organize in a more ordered structure. It is possible that 
the defects generated at the outer region of the sample 
by the sawtooth waveform diffuse deeper into the sample 
assisted by the moving lattice. As disorder is produced 
only for high ac fields (H Q > 3 Oe) that penetrate an 
important fraction of the sample, a significant portion of 
the VL will be contaminated after a certain number of 
oscillations. This description is supported by the direct 
observation in molecular dynamics simulations^ of the 
high mobility of dislocations in VL forced by oscillating 
forces. 



IV. CONCLUSIONS 

We have presented a phenomenological model that ac- 
counts for the history effects in the the ac response of 



YBa2Cu307 single crystals. The model is based on the 
effect of inhomogeneous volume-pinning in the ac sus- 
ceptibility. We obtained analytical expressions for a sim- 
ple geometry to address qualitatively the details of the 
model, and then we presented numerical results by solv- 
ing Maxwell's equations for a more realistic situation. 

We presented new experimental results that are ad- 
equate for comparison with the numerical calculations. 
The very good qualitative agreement between them indi- 
cates that we have captured the essence of what occurs 
in relation to the different dynamical responses of the VL 
inside the sample. 

The existence of a threshold sawtooth ac magnetic field 
amplitude required to lower the VL mobility supports the 
hypothesis of a ratchet-tearing of the VL produced by 
temporarily asymmetric fields as was suggested in Ref. 

Some of the differences between experimental and nu- 
merical results are probably a consequence of an intrinsic 
property of the experimental technique that undoubtedly 
perturbs the state of the VL and in a way that may de- 
pend on the initial conditions. A sinusoidal ac field in- 
creases vortex mobility and the measurement itself can 
mask the effect of a previously applied asymmetric ac 
field. Besides, J c is not restricted to take only two val- 
ues. Future experimental and theoretical work is required 
to further understand the oscillatory VL dynamics. 
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